Imaging in scattering media via the second-order correlation of light field 
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Imaging with the second-order correlation of two light fields is a method to image an object by 
two-photon interference involving a joint detection of two photons at distant space-time points. We 
demonstrate for the first time that an image with high quality can still be obtained in the scattering 
media by applying the second-order correlation of illuminating light field. The scattering effect on 
the visibility of images is analyzed both theoretically and experimentally. Potential applications and 
the methods to further improve the visibility of the images in scattering media are also discussed. 

PACS numbers: 42.50.Ar, 42.68.Mj, 42.50.Dv, 42.62. Be 



Multiple scattering has a great influence on the qual- 
ity of images. The information will be degraded and the 
images suffer reduced resolution and contrast because of 
multiple scattering, such as light propagation and imag- 
ing in the atmosphere [H, neutron imaging imaging 
and diagnosis in life and medical sciences [3| . In medical 
and clinic diagnosis, comparing with X-ray, optical pho- 
tons provide nonionizing and safe radiation for medical 
applications, and are now becoming an increasing inter- 
esting method for imaging in biological tissues. Such as 
Optical coherence tomography (OCT), Diffuse optical to- 
mography (DOT), Photoacoustic tomography (PAT)and 
so on 

E 1, S 0, S, H- Although the quality of images 
in scattering media have an enhancement in some extent 
by these techniques, there are still lots of problems which 
are difficult to be settled. Because all conventional imag- 
ing methods are based on first-order correlation of light 
field, we "see" an image only when we look at the object, 
which means the detection and imaging are unseparated 
in conventional imaging process. So when the informa- 
tion of the object is distorted by multiple scattering, and 
the information of both multiple scattering and the ob- 
ject are unknown, we can not, in principle, obtain exactly 
an image destroyed by multiple scattering. 
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With respect to the classical area of imaging, the field 
of quantum imaging aims to devise novel techniques for 
optical imaging, by exploiting the quantum nature of 
light [12J. Based on the characteristic of Bose-Einstein 
correlation of light fields and the theory of two-photon in- 
terference [lH, the imaging method by the second-order 
correlation of two light fields, for example, ghost imaging 
with entangled source or thermal light, becomes a new 
kind of imaging technique [l4| . In ghost imaging schemes, 
even if the test detector is a pointlike or bucket detec- 
tor, by measuring the second-order correlation function 
of the two light fields, we can still obtain the image of the 
object with both entangled source and thermal light by 
scanning the position of the photons which never actually 
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FIG. 1: Scheme for the second-order correlation with thermal 
light in the scattering media. 

passed through the object [H, EE EE EE EE SE SE S3 • 

This imaging method, for the first time in the history of 
imaging technology development, leads to the separation 
of detection and imaging. In this letter, when there is 
multiple scattering in the test path but no multiple scat- 
tering in the reference path, imaging in the scattering 
media with thermal light is investigated by the second- 
order correlation of two light fields. 

Fig. 1 shows the schematic of experimental setup. The 
thermal light source S produced by the method described 
in Ref. [22], SE S3, ^ rs ^ propagates through a beam 
splitter, then is divided into a test and a reference path. 
In the test path, the light goes through a thin lens of 
focal length /i, scattering media and then to the test 
detector D t . In the reference path, the light propagates 
through another thin lens of focal length / then to the 
CCD camera D r with the axial resolution 6.45x6.45/im. 

In Fig. 1, The impulse response function between the 
plane x and the plane xq is [j, Eol EH 



h(x, x Q ) = ah in (x, x Q ) + j3h aca {x, x ). 



(1) 
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where hi n (x, Xq) is the impulse response function with no 
scattering media, and h sca (x,Xo) is the impulse response 
function because of the interactions of multiple scatter- 
ing, a, [3 are the probability amplitudes of incident light 
and scattering light, respectively. 

The probability distribution function in scattering me- 
dia can be represented by point scattering function(PSF). 
Generally, there are two forms of PSF: Lorentzian-shaped 
and Gaussian-shaped distribution [2^, [26| . 



when the object plane, the thin lens f 2 and the test de- 
tector plane satisfy the thin lens equation 
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h sca (x,x ) oc J dx'P(x',x ) LA h(x,x')( L2+Z2) , (3a) 

^,,) IA .[ » (3b) 



\P(x',x ) LA \ 2 dx' = 1. 



(3c) 



z 2 + L 2 z 3 f 2 
then the impulse response function in the test path is 

h t (x t ,x 2 ) oc J dx'[ai exp | — -^a/a^j + (h J dx' 2 

xP(x', x' 2 ) LlA exp {-^4*2 }}t(x')C(x') 

xex P {^(l-£l±^)^}.(8a) 



C(x') = a 2 S(x' + m^-xt) + p 2 P(-^^x t ,x') L2 i8b) 

where t(x) is the transmission function of the object. If 
1 
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here, P(x' , xq)l a is the Gaussian-shaped point scattering 
probability amplitude from the position Xq to x' when the 
light goes through the scattering medium with effective 
thickness La- Axl a , D, w, A and n are the broaden- 
ing width, the diameter size of scattering particles, the 
concentration of suspended particles, the wavelength of 
the incident light and refractive index of the scattering 
medium, respectively, a, (3 and Axl a are all determined 
by specific experimental conditions. 

By optical coherence theory [ID, [l?], EH , the second- 
order correlation function between the two detectors is: 



AG {2 ' 2) {x r ,x t ) =< I{x r )I{x t ) > - < I{x r ) >< I(x t ) > 

dxi / dx 2 G^ 1 ' 1 \x 1 ,x 2 )h*(xi,x r )ht{x 2 ,Xt) | 2 .(4) 



and — = Z2 ^ L 2 , then the same speckle distributions can 
be obtained on the two detector planes without the scat- 
tering media and the object. Substituting Eqs. (5)-(9) 
into Eq. (4), the correlation function is 



AG {2 ' 2 \x ri x t ) oc 



x/3 1 a2}t(--j r x t) + a 1 (3 2 P{-^jX U yX r ) LlA t(yX r ) 

2 

+P 1 p 2 J dx'(x')P(x',^x r ) LlA P(-^x t ,x') L2A (10) 

If the test detector is a bucket detector, integrating xt in 
Eq. (10), then 



h h 

f XU f 

fl fl s ,/fl 



{aia 2 S(x r + x t ) + P{--jX t , -jx r ) L2A 



where <> denotes statistical average of the ensemble, 
and G tyl ' 1 \x\^x 2 ) is the first-order correlation function 
on the source plane, ht{x t ,x 2 ) is the impulse response 
function in the test path whereas h* (x r , x\ ) denotes 
phase conjugate of the impulse response function in the 
reference path. 

Suppose the light source is fully spatially incoherent, 



AG (2) (av)= J AG {2 ' 2) (x r ,x t )dx t . 



(11) 



If the test detector is a CCD camera, because the images 
on the two camera planes are inverse, in the case of x t = 
—x r , Eq. (10) can be rewritten as 



G^ 1 \x u x 2 )=I 6(x 1 -x 2 ). 



(5) 



where 1$ is a constant, and S(x) is Dirac delta function. 

Under the paraxial approximation, and when the ef- 
fective apertures of the lenses in the optical system are 
large enough, the impulse response function of the refer- 
ence system is 



AG {2a) {x r ,-x r ) oc 



otia 2 + ai(3 2 



-1 1/4 



7rA5 



2 


1/4 ^| 


ttAx 2 t 





t(j x r ) + Pi0 2 C o 



(12) 



J 71 ", 



2J7T 



h r (x r ,xi) oc exp \ jj (1 - -j)x\ - = ^-j-x r xi \ . (6) 



Co= J dx't(x')P(x',^X r ) LlA P(^X r ,x') L2A . (13) 

From Eq. (12)-(13), for i x =0 (namely, ^=0) or L 2 =0 
(namely, p 2 =0), images with high quality can always be 
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FIG. 2: (I). The axial sections of the images when the ther- 
mal light from a single slit (the silt width a=0.2mm) fixed 
before the mediums went through different mediums with the 
scattering thickness L\ + Z/2=40mm. (II). Images of the ob- 
ject obtained when the sample was fixed in different positions 
of the scattering media (averaged 15, 000 speckle frames), 
(a). Li=0mm, Z/2=40mm, Ax_L A =1.36mm; (b). Li=20mm, 
Z/2=20mm, Axl a =0.74mm; and (c). Li=40mm, Z/2=0mm, 
Aai£ A =0.01mm. (!)• the object; (2). Conventional imag- 
ing; (3)- (4). were images reconstructed by the second-order 
correlation between two paths when the test detector was a 
bucket detector, and a CCD camera with the axial resolution 
6.45x6.45/im, respectively. 



reconstructed by the second-order correlation of two light 
fields in scattering media. However, if the object is fixed 
in the middle of the scattering media, the quality of the 
image is the worst. 

In the experiments, we prepare a suspension 
liquid composed by emulsion polymerization parti- 
cles with particle diameter D=3.26/im and the so- 
lution NaCI with density p=1.19 g/cm 3 . The 
vessel used to put the suspension liquid is de- 
signed as 40mmxl0mmx20mm. We take A=650nm, 
/i=400.0mm, /=150.0mrri ) / 2 =250.0mm, z=211.0mm, 
zi=300.0mm, z 2 + L 2 =390.0mm, z 3 =243.8mm. The av- 
erage sizes of the speckles on the CCD camera plane 
are 40.6^m. When the light goes through the scat- 
tering media with thickness L\ + L2=40mm, |^| 2 «694, 
Axl a ~ 1.36mm, the scattering coefficient of the medium 
/i s ~1.64/cm. The minimum characteristic scale of the 
object ("zhong" ring) is 60/xm and the diameter of the 
ring is 1.6mm. 

Fig. 2(1) represents the point scattering functions 
when the thermal light goes through different mediums. 
The Gaussian envelop will become wider because of mul- 
tiple scattering. For conventional imaging, as shown in 
Fig. 2(H)(2), the quality of images will reduce as the 
increase of the scattering thickness L 2 . Oppositely, when 
the test detector is a bucket detector, the quality of the 
images reconstructed by the second-order correlation be- 
tween two paths will decay as the scattering thickness 
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FIG. 3: Effect of multiple scattering between the source and 
the object plane on g™(xti,xt2 = 0) and g^(x r ,xt = 0) 
for different scattering thickness L\ and 1/2=0 without the 
object, (a). Li = 10mm; (b). Li=20mm; (c). I/i=40mm. 



L\ is increased [Fig. 2(H)(3)]. However, as shown in 
Fig. 2(II)(4-c), if the test detector is a CCD camera with 
high axial resolution, we can also obtain images with high 
quality via the second-order correlation of the two light 
fields when there is only strong multiple scattering be- 
tween the object plane and the source. In addition, if 
the object is fixed in the middle of the scattering media, 
the visibility of the image will reduce, but the resolution 
doesn't decay [Fig. 2(11) (4-b)]. 

The visibility of the images reconstructed by the 
second-order correlation of two light fields can be ex- 
plained by normalized second-order correlation function 
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ff (2) (xi,a:2 = 0) = l 



AGVV(x 1 ,x 2 = 0) 
</(a?i) ><I(x 2 =0) >' 



(14) 



where {x\ , x 2 = 0) denotes the second-order cross cor- 
relation degree between two spatial positions x 2 =0 and 
x\. For the thermal light field, g^ 2 \xt\,xt 2 — 0) reveals 
the fluctuation of the light field, whereas g( 2 '(x r , Xt — 0) 
describes the cross correlation between two paths. As 
shown in Fig. 3, both the maximum values of the 
second-order correlation degree g( 2 \x t i = 0,x t2 = 0) 
and g( 2 \x r = 0,x t — 0) decrease sharply as the scat- 
tering thickness L\ is increased, so the visibility of im- 
ages will reduce (see Fig. 2(H)(3)). However, even 
g( 2 \x r = 0, Xt = 0) is lower than 1.05, an image with 
high quality can still be reconstructed when both the de- 
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tectors in two paths are CCD cameras with high axial 
resolution (see Fig. 2(H)(4)). 

Because multiple scattering between the source and the 
object plane destroys the cross-correlation between the 
two paths, thus the visibility of the image and g^ (x r — 
0,xt = 0) will be degraded as the increase of multiple 
scattering. Higher g( 2 > (x r = 0, x t = 0) and larger trans- 
verse coherent width on the object plane can both en- 
hance the visibility of the image in scattering media re- 
constructed by the second-order correlation between the 
test and reference paths [13, HI] . 

In addition, by the results described by Eq. (12) and 
represented in Figs. 2-3, when the both detectors in two 
paths are CCD cameras with high axial resolution, the 
quality of images obtained by the second-order correla- 
tion of two light fields is also determined by the quality 
of images registered by the CCD camera D t . Thus, al- 
most all existing conventional imaging schemes in scat- 
tering media can be applied in the test path to further 
improve the quality of the correlation imaging system. 
When Li « L 2 , the quality of the image reconstructed 
by the second-order correlation of the two light fields is 
mainly determined by the mechanism of "ghost" imag- 
ing. Because of the separation of detection and imag- 
ing for "ghost" imaging, the test detector is only used 
to collect the radiation of the thermal light propagating 
through the object, while imaging system is located in 
the reference path. Because there is no multiple scatter- 
ing in the reference path, the quality of the image isn't 



influenced. If L\ >> L 2 , the quality of the reconstructed 
image mainly depends on the conventional imaging sys- 
tem in the test path. Because the scattering between the 
object and the CCD camera D t is week, thus the im- 
age with good quality can also be obtained. Generally 
speaking, by the second-order correlation of the two light 
fields, we can always obtain an image with much better 
quality than the image achieved only with conventional 
first-order correlation optical imaging methods in scat- 
tering media. Entangled source and other nonclassical 
light source with higher g^{x t \ — 0,x t2 = 0) may be 
applied to further improve the quality of the image ob- 
tained by the second-order correlation imaging system in 
scattering media [1^, [3(3, HU . 

In conclusion, imaging via the second-order correlation 
of two light fields provides a brand-new route for imag- 
ing in scattering media. We demonstrate for the first 
time that we can always obtain an image with much bet- 
ter quality reconstructed by the second-order correlation 
imaging method than the image achieved by conventional 
first-order correlation imaging methods in scattering me- 
dia. This will be very useful to imaging, test and diag- 
nosis of biological tissues with infrared and near infrared 
light. 
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